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Abstract 

We investigate the implications of the long-rang self-interaction on both the self-capture and 
the annihilation of the self-interacting dark matter (SIDM) trapped in the Sun. Our discussion is 
based on a specific SIDM model in which DM particles self-interact via a light scalar mediator, 
or Yukawa potential, in the context of quantum mechanics. Within this framework, we calculate 
the self-capture rate across a broad region of parameter space. While the self-capture rate can 
be obtained separately in the Born regime with perturbative method, and in the classical limits 
with the Rutherford formula, our calculation covers the gap between in a non-perturbative fashion. 
Besides, the phenomenology of both the Sonnnerfeld-enhanced s- and p -wave annihilation of the 
solar SIDM is also involved in our discussion. Moreover, by combining the analysis of the Super- 
Kamiokande (SK) data and the observed DM relic density, we constrain the nuclear capture rate of 
the DM particles in the presence of the dark Yukawa potential. The consequence of the long-range 
dark force on probing the solar SIDM turns out to be significant if the force-carrier is much lighter 
than the DM particle, and a quantitative analysis is provided. 
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I. INTRODUCTION 


Although the existence of dark matter (DM) have been well established by the observa¬ 
tions from the galactic scale up to the cosmological scale, the particle nature of DM still 
remains unclear. While the collisionless Weakly Interacting Massive Particle (WIMP) lias 
been the focus of the current DM study, the possibility of the self-interacting dark matter 
(SIDM) has also drawn increasing attention in both particle physics and astrophysics (see 
e.g. |T| 0| ) ■ Long-range self-interactions between DM particles induced by the exchange of 
light mediator particles can lead to the Sommerfeld enhancement of DM annihilation cross 
sections fz-12], which can be used to explain the observed cosmic-ray positron excesses 13- 


16] by providing a source of the boost factor of (P(10 2 — 10 3 ) (for a recent global analysis on 


the boost factor, see e.g. 17h23|). It is under active investigation whether an excess of an¬ 


tiproton also exists in the AMS-02 p/p data 24-131) and a boost factor is also required in the 
antiproton channel. Phenomenological consequences of the light mediator exchange between 


DM particles in astrophysics and cosmology have been discussed in a series of studies 10 


0.0 

49 

shape 

35, 


which include its implications on relic density (HI, 0, B, S3 , CMB IHS] , halo 
42, (50], and small scale structures [32], ]4g, |51|, etc. 


Another interesting phenomenology associated with the DM self-interaction is its self¬ 
capture within the Sun. In addition to solar elements, DM particles that have already 
been captured within the Sun may also contribute to the capture of the passing halo ones. 
Consequently, more DM particles will build up within the Sun and an enhancement of 
the high-energy neutrino signals from the DM annihilation might be possible. Under the 
assumption of a constant self-scattering cross section, the authors of Ref. 52-5^] utilize 
neutrino detectors to probe the solar SIDM particles and explore for the viable parameter 
space. It is found that a significant enhancement of the solar neutrino flux requires a 
disparity between two separate input parameters, i.e., a large DM self-capture rate and a 
small annihilation cross section. In this paper, we will extend this case to the scenario where 
the DM particles interact with each other via a light mediator (also in Ref. [551,156]*). In such 
scenario, this force-carrier will play a role in both the self-capture and annihilation within 
the Sun, as well as in thermal freeze-out that accounts for the observed DM density. As a 
result, the self-capture and annihilation rate become inter-related to each other, and more 
than that, the disparity between them may be achieved in a natural way. To investigate the 
relevant implications from an optimistic view, we examine a concrete SIDM model in which 
the Majorana DM particles y interact with each other through the exchange of a light scalar 
mediator 0, so that the “irreducible” annihilation process yy ~t </></> proceeds through the p- 
wave channel. The p-wave annihilation is of particular interest because contrary to its s-wave 
counterpart that is subject to the Sommerfeld enhancement, it suffers a velocity-suppression 
not only in the Sun but at the kinetic decoupling stage, and hence results in a larger relic 


density 


12 


As a consequence, a larger self-interaction coupling a x is required accordingly 


for the p-wave annihilation to give the correct relic density, which is also a prerequisite for 


In these two papers, the self-capture rate is calculated either in the Coulomb limit [55] or with the pertur¬ 


bative method 


56]. 
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111. 135l . 42] in which authors managed to discuss the DM self- 


a large self-capture rate. 

Following the literature 
scattering beyond the perturbative level, we similarly give a quantitative description of the 
self-capture rate for the Majorana DM particles. However, instead of adopting the transfer 
cross section as a proxy quantity to describe the self-scattering, we calculate the self-capture 
rate with long-range interaction in a straightforward manner, considering that the kinetic 
requirement for self-capture automatically regulates the forward-scattering divergence and 
hence makes the self-capture rate a well-defined quantity. 

In addition, in order to explore the phenomenology of the solar SIDM from a compre¬ 
hensive perspective, our discussion includes the Sommerfeld-enhanced annihilation of the 
solar SIDM as well. We illustrate that while the DM particles predominantly annihilate to 
light mediators through the p -wave channel, a small but non-vanishing quota of annihilation 
proceeding through the s-wave channel in the thermal freeze-out may become dominant 
over its p -wave counterpart within the Sun, and produce signals at the terrestrial neutrino 
telescopes. For such case we use the neutrino telescope Super-Kamiokande 57] to find out 
whether the Sommerfeld effect can help constrain the nuclear capture rate C rt in a nontrivial 
manner. These constraints can be translated into upper limits on the DM-nucleon cross 
section for other SIDM models. It is found that in the approach of the Coulomb limit of the 
long-range self-interaction, the constraints from the Super-Kamiokande on nuclear capture 
rate can be significantly tightened. 

This paper is organized as follows. In Sec. [TT] we take an introduction to the model 
setup for our study. In Sec. IIIII we introduce the related background about the solar DM, 
including its self-capture and evolution within the Sun. We discuss the implications of 
the light mediator on the self-capture and annihilation of the solar DM particles, as well 
as on relic density in Sec IIV1 We further investigate its impact on the neutrino telescope 
Super-Kamiokande in Sec. El and conclusions are drawn in Sec. ED 


II. MODEL SETUP 


We are interested in the the dark sector that gives rise to the Yukawa potential of the 
DM self-interaction. As a familiar prototype, Yukawa potential between DM particles pro¬ 
vides us with a representative description and a concrete example of the long-range self¬ 
interaction lo], EsiH]. The self-interacting Yukawa potential* is written as 


V (r) = ±^ e ~ m * r : 
r 

which describes in non-relativistic limit the following interactions of the dark sector 


( 1 ) 


42 


{ QxXl^X^ti vector mediator 
QxXX^ scalar mediator 


( 2 ) 


Strictly speaking, Yukawa potential narrowly refers to the attractive interaction in Eq. 0. However, for 
the purpose of convenience, here we use it to refer to both the attractive and repulsive interactions. 
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Figure 1: The tree-level diagrams for elastic scattering process XX XX 


where m$ is the mass of the mediator particle 0, and the dark fine-structure constant 
a x is connected to the coupling g x through a x = g x / ( 4 it). For a vector force-carrier, 
the interaction between DM could be attractive (xx scattering) or repulsive (xx. or XX 
scattering), while the force can only be attractive (—) between DM for the scalar mediator 
scenario. Since in this paper we are interested in spin-1/2 identical DM, i.e., Majorana DM 
particles, the interaction between DM and the vector mediator in Eq. (J2J) is absent, and an 
extra s-channel process should be included in the calculation of relevant scattering amplitude 
for a scalar mediator </>, which is shown in Fig. CD While the second and the third diagrams 
in Fig. Q] correspond to the long-range Yukawa interaction in non-relativistic limit, the first 
s-channel process can be described non-relativistically with a contact interaction, which is 
subject to the Sommerfeld enhancement that will be reviewed in detail below. However, due 
to its highly suppressed cross section o xx ~ a x v^/m x , contribution of the s-channel process 
can be safely ignored compared with that of the Yukawa potential, even if the Sommerfeld 
enhancement is involved. So only the Yukawa potential is relevant for the DM scattering. 

Although we want to explore phenomenologically the impact of the DM long-range self¬ 
interaction as far as possible, it is unrealistic for our discussion to be entirely model- 
independent. At the very least, the phenomenological consequence of the irreducible an¬ 
nihilation process XX 00 has to be taken into consideration (see Fig. [2]). So, we study 
a model that not only gives rise to signihcant Sommerfeld effects within the Sun, but also 
cover as wide a range of phenomenology as possible, namely, the self-capture and both the 
s- and p-wave annihilation of the solar SIDM. 

To facilitate the novelty brought by the long range self-interaction, we assume that the 
irreducible p-wave process (shown in Fig. [2j) dominates the annihilation during the ther¬ 
mal freeze-out. The reason for this assumption is two-fold: 1) As opposed to the s-wave- 
dominated case, the p-wave-dominated annihilation avoids severe constraint on the coupling 
a x imposed by the DM relic density, so that the Sommerfeld effect remains significant 12, [58|. 


2) The p-wave-dominated annihilation within the Sun is subject to such signihcant velocity- 
suppression that the total annihilation rate will be enhanced remarkably (see Eq. (J9]) 
and (TT8j) in the following section), which may further favor the indirect detection of the 
solar DM. 

On the other hand, given the dark sector of the Majorana DM x and the scalar mediator 
0, it is natural to assume that they couple to the SM sector through the Higgs portal 59 65|. 
In addition, extra four-fermion interactions are introduced to account for the DM s-wave 
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Figure 2: The diagrams for the minimal annihilation process XX 00- 


annihilation into the SM particles. Integrating all these considerations, we consider the 
following interaction as an example: 


£ => y xx4> - y f sin 6 + y^ N °N + xi 5 xff , 

/ / 


( 3 ) 


where sin 0 describes the mixing between 0 and Higgs boson, and yj — irtf/v ew is the 
Yukawa coupling of the SM fermion /, with the SM fermion mass m/ and the Higgs vacuum 
expectation value vew ~ 246 GeV. Gf is the coupling strength between the DM and the SM 
fermion /. The coupling of 0 and right-handed sterile neutrino N is introduced in Ref. {66] 
to save the light mediator 0 from the Big Bang Nucleosynthesis (BBN) constraints. To 
avoid the overclosure problem, overproduction of entropy and changes to the light element 
abundances during the BBN, 0 is required to decay before the BBN (t Is) jfi3], which 


corresponds to the bound 66 


Un <1 6 x 10 


-12 


100MeVy /2 


m ^ 


( 4 ) 


III. DARK MATTER IN THE SUN 
A. Dark matter in the Sun 

The accumulation of the DM in the sun is initiated by DM-nuclei collisions that cost the 
passing-by Galactic DM particles sufficient kinetic energy to escape from the gravitational 
pull of the Sun. In this process, subsequent collisions tend to further lower the kinetic 
energy of the DM rather than eject them back into the deep space as long as the DM mass 
is heavier than a few GeV, hence the number of the trapped DM grows until their capture 
and annihilation reach an equilibrium. If the DM self-interaction is taken into account, the 
evolution of the solar DM number N x is depicted by the following equation 

^ = C n + (G s - E n ) N x - (C a + E s ) Nl (5) 

where C n is the DM capture rate by (scattering off) nuclei in the Sun, C a is twice the 
annihilation rate of a pair of DM particles, E n and E s represent respectively the evaporation 


5 








rate by nuclei and by the captured DM particles, and C s is the capture rate by the trapped 
DM particles, a key physical quantity of our concern in this paper. It should be noted that 
in the range of the DM mass of our interest (m x > 10 GeV) it is justified for us to omit 
the the evaporation rate E n in Eq. (|5j) , since it is irrelevant for a DM mass larger than 
5 GeV 68 h 70|.* Besides, we also omit the the self-evaporation rate E s because C a E s , 
which will be demonstrated in detail in Appendix [B] Thus Eq. ((5]) is reduced to 


dN x 

dt 


= C n + C S N X - C a N% 


which has an analytic solution 


N x = 


C n tanh (t ■ £) 
f - (Cs/2) tanh (t ■ £)’ 


( 6 ) 


(7) 


with 


£ = y/Cn-C> + Q/4. ( 8 ) 

If tanh (t & ■ £) ~ 1 in Eq. ([7]) with solar age t & = 4.5 x 10 9 yr , the capture and annihilation 
of the DM reach an equilibrium today. Then the annihilation rate r A for a steady number 
of DM within the Sun can be written as 

r A = \c a Nl 


= ,rG 



“ l\ c * + ^k + \ k + uk' Cs 


(9) 


which is reduced to the trivial case where the dark matter self-interact ion is absent, i.e., 
when r A = C n /2 if C^/(C n C a ) -C 1 is satisfied. 

B. Self-capture and annihilation of the self-interacting dark matter 


The formula for the _self-capture rate can be obtained after a minor modification to the 
procedures in Ref. 69, 72] as the following: 

C S N X = — f n x (r) d 3 r [ a sc (w) —f (u) d 3 u, (10) 


m 


X J Sun 


U 


where p x = 0.3 GeV is the DM local density in our solar neighborhood, m x manifestly the 
mass of the DM. / (u) is the velocity distribution of DM in the Sun’s rest frame, which is 
expressed as 


/( u ) = 


( u+v q) 
e 3 

IV(Uesc) ’ 


( 11 ) 


For smaller DM masses, the detection of the solar DM evaporation is discussed in Ref. [71]. 
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where |v 0 | = 220 km • s' 1 is the velocity of the Sun, u being the DM velocity at infinity where 
the solar gravitational effects is negligible, t>o = 220 km • s -1 the DM velocity dispersion, 
and N(v esc ) is the normalization constant dependent on the Galactic escape velocity v esc = 
544 km • s _1 . w = s/vl sc (r) + u 2 is the velocity of the incident DM particle, with v esc ( r ) 
the solar escape velocity at radius r. The DM distribution within the Sun can be described 
with a characteristic temperature T x as follows* 


n x (r) oc exp 




GM (r') 

r t 2 



( 12 ) 


where G is the Newton’s constant and M (r) represents the solar mass contained within 
radius r. T x is used to be determined self-consistently through the following equation: 


T x — T q (R x ), 


(13) 


in which the right hand side represents the solar temperature at the mean value of DM 
radius 


f Sun n x ( r)rd 3 r 
X Jsun n x (r) d ‘R ' 


(14) 


We use the Standard Sun Model (SSM) GS98 [73] to obtain a fitting function of R x as 
follows _ 


R x ~ 0.012 x 


'100 GeV 


771 , 


R& , 


(15) 


with solar radius R & . The effective cross section for self-capture, cr sc ( w ), is dependent on the 
velocity of the incoming DM particle w, if we approximate the target DM trapped within the 
Sun to be at rest compared with the accelerated Galactic ones. On the other hand, for the 
mass range concerned in this paper (m x > 10 GeV), we further assume that the captured 
DM particles are located at the center of the Sun (see Eq. (fT5l) h thus the self-capture rate 
can be simplified as 



where w = = \A ; esc (0) + u2 ■ 

However, here we note that instead of growing infinitely with the accretion of the captured 
DM particles, C S N X saturates to a maximum value set by the geometric limit of the resident 
DM, above which C S N X is replaced by 


(C S N : 


XJ max 


Px 

m x 

Px 


m, 


tcR , 


7r r : 


vie ^ 


u 

vl + u 2 

U 


f (u) d 3 u 


/( u) d 3 U 


~ 4.67 x 10 24 x 


lOOGeVy 


-i 


m. 


To avoid distraction we postpond a discussion on the solar SIDM distribution to Appendix [A] 


( 17 ) 
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On the other hand, the annihilation rate of the solar DM can be expressed as 


C a 


(^anriD-e^Q 


IsunK ( r ) ^ 

C/sun n x ( r ) d 3 rj 2 ’ 


(^"anrTrel) q 

VeS 


( 18 ) 


where (cr ann n re i) e is the thermal average of the DM annihilation cross section times relative 
velocity n re i within the Sun, and the effective volume occupied by the trapped DM particles 
can also be described with a fitting function as follows: 


D eff = 6.9 x 10 27 


/100 GeV\ 3/2 

V m x ) 


cm 3 . 


(19) 


IV. IMPLICATIONS OF THE LONG-RANGE SELF-INTERACTING DM 

In this section, we study the implications of the long-range self-interact ion on the self¬ 
capture and annihilation within the Sun. Due to the multiple exchanges of the light force- 
carrier between DM particles, the relevant calculations have to be carried out beyond the 
perturbative approach. The quantum mechanical description of the two-body scattering is 
proved to be an effective framework to do so, which is connected to the self-capture and 
annihilation in a manner that the asymptotic wave-function at infinity describes the final 
state of the scattering, and the wave-function at origin partly determines the probability of 
the annihilation. Now we delve into the details. 


A. Calculation of the self-capture rete 


We calculate the self-capture cross section cr sc ( w ) in Eq. (TTTTh in the center-of-mass (CM) 
frame, in which the self-capture differential cross section for a pair of identical DM can be 
expressed as (after averaging over the initial spins and summing over the final spins ) 


da 


SC 


dVt 


1 

2 


E< 2S+1 ) 

. s 



/((?) + (-1) S -/(7T - 



( 20 ) 


where S represents the total spin of the DM-pair, the factor 1/2 accounting for the symmetry 
of the out-going DM particles, and the elastic scattering amplitude / ( 9 ) is related to the 
wave function in position space via 


< x I^p> = 

large r 


l 

(2tt) 3/2 

1 

(2tt) 3/2 


Y, i e (2f + 1) At (r) Pi (cos 9) 

i 

pipr " 

e ipx + /(0)— , 

r 


( 21 ) 
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m =10 GeV 



Figure 3: Shown is the self-capture rate C s dependent on parameter z = m x /m ( j > and a x for the 
Majorana DM with mass rn x = 10 GeV. 


where the radial partial-wave amplitude Ap (r) has the following asymptotic solution for a 
large r: 

A e (r) = e lSl [cos Sp jp ( pr ) - sin 5p n £ ( pr )], (22) 

with je(np) as the spherical Bessel (Neumann) function, and Sp the phase shift for a partial 
wave t . Thus / (9) can be further expressed as 


1 


+oo 


f (0) — - y (2t + 1) e lSe sin SpPp (cos 9 ), 
^ e=o 


(23) 


with p = m x w/ 2. As for the scattering angle 9, to ensure a net capture it is required to 
satisfy 

— c < cos 9 < c, (24) 


with 


u 


c = 


2 y j 


+ u 2 


(25) 


where fj = 18.5 km • s^ 1 is the solar escape velocity at the radius of Jupiter’s orbit. As a 


conservative estimate |74|, uj is introduced not only to take into account the gravitational 
effects of the Jupiter from a more realistic three-body perspective, but to regulate the 
divergence encountered in the calculation of C s for the long-range self-interaction scenario. 

We follow the same procedure as proposed in Ref. 42] to numerically calculate the phase 
shift Sp and then the self-capture rate C s for specific parameters m x , m and coupling 
strength ot x . The value of C s is scanned in Fig. [3] over the the parameter space of the 
mass ratio z = m x /m^ and the self-interaction strength a x , for the Majorana DM with 
mass m x = 10 GeV. While the Born regime occupies the lower left corner of the figure, the 
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Sommerfeld effect becomes increasingly significant in the Coulomb limit (z —>■ oo), where a 
large partial wave number (£ max > 1000) is required in calculation to ensure the convergence 
of C s . It is worth mentioning that given z and a x hxed, C s will be suppressed for larger DM 
masses through a m“ 3 dependence, which is obvious from Eq. (1T6l 1201 and [23]) . 


B. Sommerfeld enhancement of the annihilation rate 


Apart from the self-capture, the light force-carrier brings about significant consequences 
on the annihilation of DM as well, which have been widely studied in determining the thermal 
relic density and in calculating signals in indirect DM searches. To gauge the Sommerfeld 
effect in annihilation we invoke the Sommerfeld factor as follows 


75, 76 


5, = 


non—perturbative partial wave cross section 
tree — level partial wave cross section 


M l 


ladder 


M 


(26) 


where ^ stands for the £-th partial wave. If one skips somewhat formal theoretical discussions, 
the implication of Sommerfeld enhancement on annihilation is encoded in the following 
heuristic relation between the non-perturbative ATf adder and “bare” amplitude Mq at tree- 
level 


75, 76 


M 


ladder 




(27) 


where (j) (k) = (k ) + (—l) 5 k ) /2 is the familiar scattering wave function for 
identical DM-pair in the reduced system, p = m x v re \/2 being the initial momentum in the 
CM frame, with relative velocity v re i. It is noted that for identical fermion pairs L + S must 
be even, thus we have </> (k) = (k |-0p ). After Fourier transformation and by repeated use 
of addition theorem the Sommerfeld factor for partial wave i can be written in terms of the 
wave function in position space 


5, = 


(2£+l )\\d l A t (r) 


p f t\ 


dr 1 


r =0 


(28) 


which can be calculated numerically along with 8g. Therefore, after taking into consideration 
the Sommerfeld enhancement, we express the annihilation rate for the trapped DM particles 
as follows 


C a = 


(^ann^rel),- 


0 




= / ((Tann |v 2 - Vj|) tree <S / x (v^ f x (v 2 ) 

= y~s I ^ crannWrel ^ tree ^ 


(29) 


where S is the Sommerfeld factor dependent on u re i, the footnote “tree” indicating that 
the cross section is calculated at the tree-level, and the Maxwellian velocity distribution is 
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Figure 4: Shown are the thermally averaged cross section (a-, 
and p -wave annihilation within the Sun, respectively. See text for details. 


ri re i) e . with a x = 0.01, for the s-wave 


written as 


fx ( v ) = 


e -v 2 /ul 
2 ) 3 / 2 : 


(30) 


7 TU 


with u 0 = sj2T x /m x . Note that in Eq. (l30l) we neglect the cut-off at the solar escape 
velocity, which is equivalent to the approximation v esc —> oo. 

Since we are only interested in the case where other p-wave processes are suppressed by 
couplings y /sin 6 or y n in Eq. (l3lh the total thermally averaged cross section can be written 
as 




i l Vel) ((^^reOtree *^l) + ( <T // Vrel ) tree ^ 

= ((<h^ re i) tree 5i> + (p (crn rel ) 0 S 0 ), 


where included are the thermal cross sections for the irreducible p -wave process 

{ a H V rel) tree 


37ra 2 u 2 el 


8m?, 

A 


and its s-wave counterparts arising from the four-fermion interaction 

q2 / ,™2 \ 3/2 




(31) 


(32) 


(33) 


with the color factor c/ = 3 (1) for quarks (leptons). <Si (So) is the p -wave (s-wave) Som- 
merfeld factor. The s-wave component is parametrized as a small factor p times the typical 
averaged thermal cross section (<7U re i) 0 = 3 x 10 -26 cm 3 • s -1 in the second line of Eq. (J3TJ). 
To illustrate the resonant behaviors of the annihilation, in Fig. [4] we present the solar ther¬ 
mally averaged cross section ((cry re i) 0 Sq) q and ((cr^n re i) tree <Si) Q as the functions of the DM 
mass m x and z = m x jm ^, respectively, with coupling a x =0.01. The DM mass parameter 
m x actually reflects the dependence of the DM velocity within the Sun through the relation 
u 0 = yfiTJ m x in Eq. 
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Figure 5: The coupling strength a x required to give the correct relic abundance as a function of 
DM mass m x for z = 10 (black dashed) and 10 3 (red solid), under the condition g = 10 -2 . See 
text for details. 

C. Constraints from the relic density 


In Sec JIV Al we calculate the self-capture rate C s with free parameters a x and z = m x /m 
However, the coupling strength a x is not an independent parameter as it can be determined 
by the observational DM relic density with specific m x and m^. For the given m x , m$ 
and a x the evolution of the relic density of DM is described by the following Boltzmann 


equation 77], 178 


dn x 

dt 


+ 3 Hn x - - (cr ann u re i) (n 2 - n^ 2 ) , 


(34) 


where n x and n® q are the number density of the DM and its counterpart in equilibrium with 
the background particles, respectively, H the Hubble constant and (caimD-ei) is the thermally 
averaged cross section introduced in Eq. 0311) . In practice by changing variables such that 
y = (Y‘ q = 'n® q /s) with s the entropy density and x = m x /T , the above Boltzmann 
equation can be rewritten as follows |77l . |78| 


dY 

dx 



7T (g*s/\JH*) 

mpm x - --- 

X 2 


(dannDel) ( Y 2 - F eq2 ) , 


(35) 


where g* s and g * are the effective degrees of freedom for entropy and energy density, respec¬ 


tively, and the Planck mass rrip = 1.22 x 10 19 GeV. We follow Ref. 58, 79] to numerically 
solve Eq. (j35l) so as to give the coupling a x determined by the relic density. In Fig. [5] shown 
is the value of a x as a function of m x and 2 = m x /m with a specific value rj = 10” 2 intro¬ 
duced in Eq. (13TT) . One can see from Fig. [5] that owing to the high velocity at the chemical 
decoupling stage, a x is insensitive to the Sommerfeld effect. Here we note that the effect 
of the kinetic decoupling is not included in our calculation, because it is insignificant for 
the case of p- wave-dominated annihilation when the temperature of kinetic decoupling Thd 
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Figure 6: The solar thermal cross section (cranntyei}^ determined by the relic density for the s-wave 
and p -wave annihilation, respectively. See text for details. 


is much less than that of the chemical decoupling Tf, even if the Sommerfeld enhancement 
is taken into consideration [42l. l58l. 8()|. This requirement is in coincidence with the onset of 
the Sommerfeld effect: a large z = m x /m$ 1 guarantees that Tf ~ m x / 25 ~ T^ d , 

where the kinetic decoupling temperature for the y-0 interaction T^ d provides a maximal 


value for 12|. Given 7] = 0.01, in Fig.[6]we use the relic density to calculate the thermally 


averaged cross section for the s-wave ((cry re i) 0 <S 0 ) Q and p -wave annihilation i) tree «Si) 

within the Sun, respectively. We note that smaller r) will not bring remarkable changes to 
Fig. [6] because the relic density turns out to remain constant as the s-wave component gets 
smaller (r) < 0.01) fhsl j. 


V. CONSEQUENCES ON THE NEUTRINO TELESCOPE 


The signal of the solar DM is associated with the observation of primary and secondary 
high energy neutrinos, the annihilation products of the DM accumulated in the center of 
the Sun. Among other goals, many terrestrial neutrino detection projects such as Ice- 
Cube [j81j, Super-Kamiokande [82j, Baikal Neutrino Project [ssj and ANTARES are 
dedicated to such observation, and next generation telescopes including PINGU |85| and 
Hyper-Kamiokande 86] are also in proposal. In general the neutrino differential flux at the 
detector location can be schematically described as 


d&v 

dE„ 


dN,, 


47t dE u 1 


(36) 


where is the annihilation rate introduced in Eq. (I51h d Q the Sun-Earth distance, and 
dN u /dE u denotes the energy spectrum of neutrino u e:AtiT and P eiAt)T per DM annihilation 
event, after folding the effects of annihilation branch ratio, hadronization of quarks and 
neutrino oscillation, which are usually simulated with DarkSUSY 87] and WimpSim 83]. For 
the purpose of illustration the neutrino flux is connected to the readout at the detector 
terminal with the following relation: 
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(37) 




dt 




d<S> v 

dE,. 


A u (E v ) dE u , 


J 0 «/ J -^th 

where A0 is the number of the muon events that arise from the charge and neutral current 
interaction between the incident neutrino and medium surrounding the detector, the effects 
of which is encoded in the effective area A v , E t h the threshold energy of the detector, 
Ai? = 27t(1 — cosT) being the solid angle surrounding the Sun with angular resolution T, 
and t exp is the exposure time. If no excess of muon over the background is observed, an 
upper limit on the annihilation rate IA can be obtained. 


Now we use the recent release of the Super-Kamiokande (SK) analysis |57| to constrain 
the strength of the DM-nucleon interaction in the presence of the Sommerfeld effect, in terms 
of the upper limits on the nuclear capture rate C n that depends on the specific DM-nucleon 
interaction model. One can translate these constraints back on the couplings of the DM- 
nucleon interaction with the nuclear form factors given by Ref. 89, 9oj, once the interaction 
is specified. For the purpose of simple illustration we assume a leptophilic four-fermion 
interaction in Eq. (J3j), so that the incoming DM particle with velocity w interacts with the 
target nucleus through the Higgs portal. The relevant differential cross section is expressed 


as 


d° X N c2 ■ 2 n 

= «x/ W sm 9 


m p \ 2 A 2 F^(q 2 ) 1 


Tew / (q 2 + 


m 


2^W 2 


(38) 


where q is the transferred momentum of the incident DM particle, and A is the atomic 
number of the target nucleus. The nucleon mass is approximated as the protonic one m p , and 
the coefficient fA = 0.351 arises from the scalar bilinear for the underlying quarks. F 2 n (Q 2 ) = 
[3 ji ( qRi ) / (q R\)] 2 e~ q2s2 is the nuclear form factor, where ji(x) = sin (x) jx 2 — cos (a:) jx 
the Bessel spherical function of the first kind, and the parameter R\ = \JR^ — 5 s 2 with 


is 


R 0 ^ 1.23H 1//3 fm, and s ^2 1 fm 911. 

Current DM direct detections provide a stringent upper bound on the mixing angle such 


that sin 9 < 10 8 [66|. Mixing angle at such scale is 2 to 3 orders of magnitude smaller 
for the light force-carrier 0 to decay before f ~ Is, jeopardizing the success of the BBN, 


so 0 is assumed to decay quickly to the sterile neutrino N instead [66|, as explained below 
Eq. (J3]). As a consequence we regard the s-wave annihilation as the only source relevant 
for the DM indirect detection. In Fig. [7] we present the upper limits at 90% confidence 
level (C. L) on C n imposed by SK with various z = m x jm^ = 10 2 , 10 3 , andlO 4 , and 7] = 
10~ 2 , 10~ 4 , 10 -5 andlO -6 , respectively. As is discussed similarly in Ref. [57;], two extreme 
benchmark scenarios in which DM exclusively annihilate via hard (t + t“) and soft (bb) s- 
wave channel are taken into consideration over a DM mass range from 10 to 100 GeV. The 
number of the trapped solar DM particles N x , or equivalent ly. the nuclear capture rate C n 
is constrained by the upper limit on the annihilation rate T up with the following relation: 


Eu p A 2 C a N 2 B¥L s , 


(39) 


where the branch ratio of the s-wave annihilation which is responsible for the neutrino signal 
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Figure 7: SK 90% C. L. upper limits on nuclear capture rate C n over a range of DM masses for 
rj = 10 —2 , 10 -4 , 10~ 5 andl0~ 6 , respectively. Solid (dashed) lines indicate hard (soft) annihilation 
channels. Shown in the dash-dot line is the nuclear capture rate calculated with parameters sin 6 = 
10 -8 and z = 10 4 . See text for details. 


can be expressed as 


BR, = (ri(av, el ) 0 S 0 ) / + (i)(<Tti r6l ) 0 5 0 ) 


(40) 


The constraints on C Q are affected mainly by three factors: T up , z and rj. As a rough 
estimate, Eq. (J9]) can be used to explain the features of the constraints in Fig. 0 when the 
Sommerfeld effect is insignificant {z = 10 2 , 10 3 ) so that T up C%/ (2 C a ), upper limits are 
sensitive to rj through the dependence on BR S ; the sharp increases of the sensitivity on C n 
in Fig. 0 imply that r up ~ O (C% BR S / (2 C a )) when the Sommerfeld effect turns remarkable 
(z = 10 4 ), which is of our interest in this study. For reference purpose in Fig. 0 we also 
present the nuclear capture rate (in the black dash-dot line) derived from the DM-nucleus 
cross section in Eq. ()38l) . for which the values of the 0-Higgs mixing angle sin 9 = 1CT 8 
and z = 10 4 are adopted. Parameters at such scales are comparable to the sensitivity of 
direct detection discussed in Ref. [66]. Although the direct detection approaches are proved 
to be efficient in exploring the spin-independent nucleon-DM interaction, the solar neutrino 
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Figure 8: For the given 2 = 10 3 and z = 10, branch ratio of the s-wave channel as the function of 
m x in the Sun. 

detection may impose even more stringent constraints on that as the relevant parameter 
space is significantly squeezed by the Sommerfeld effect. Such phenomenon occurs at z = 10 4 
for the present SK sensitivity, as is shown in Fig. [71 

It is also noted that when the s-wave component is sufficiently small that its cross section 
drops below the p-wave counterpart, constraints will turn loose because that’s when BR S 
begins to be suppressed by a small p. To get some sense, in Fig. [8]we plot BR, as a function 
of DM mass m x for z = 10 3 and z = 10 4 . It is evident that the dips in BR S in the left panel 
are attributed to the p-wave resonance encountered along the cross section at z = 10 3 in 
Fig. [6j and hence are responsible for the corresponding bumps in sensitivity on C n in Fig. [7] 
for z = 10 3 . 

VI. CONCLUSIONS 

In this work we investigate an example of the solar DM that self-interact via light force 
carriers, specifically, Majorana DM particles scattering with long-range self-interaction and 
annihilating dominantly through p-wave channel. To optimize and cover the relevant phe¬ 
nomenology of the solar SIDM as comprehensively as possible, we study a scenario where a 
small but non-vanishing quota of DM annihilation proceeded through the s-wave channel in 
the thermal decoupling epoch, so that the self-capture, both s- and p-wave annihilation are 
all involved in our consideration. 

Instead of evaluating relevant DM self-capture and annihilation processes at tree level, we 
calculate them in a non-perturbative way by solving the Schrodinger equation with Yukawa 
potential. We find that the Sommerfeld effect can be significant for DM self-capture if 
coupling a x and mass ratio z are regarded as free parameters. Besides, the impact of the 
light mediator on the thermal freeze-out that is assumed to produce the the observed DM 
relic density is also taken into consideration, and is used to determine the self-interaction 
coupling a x for the given masses of DM particle and mediator. 
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Moreover we also explore the consequence of the Sommerfeld effect on terrestrial neutrino 
telescopes which is utilized to probe signals from the solar DM. Since the light mediators 
produced from the p-wave DM annihilation do not produce any observable signal at the 
terrestrial neutrino telescope, the neutrino signals originate exclusively from the s-wave 
component, which may become close to or even dominant over its p -wave counterpart within 
the Sun, due to the Sommerfeld enhancement. Therefore we combine the sensitivity of SK 
and the observed DM relic density to constrain the nuclear capture rate, or equivalently, 
the strength of DM-nucleon interaction. Due to SK’s high sensitivity, the constraint on the 
DM-nucleon interaction can be remarkably tightened in the low DM mass region, where the 
Sommerfeld effect on the DM self-capture rate has not yet been outstripped by the m^ 3 
suppression. 

It is also noted that while our discussion is based on a generic class of models in which 
the fermion DM particles self-interact via a light scalar mediator, the consequence of this 
force-carrier in astrophysics and cosmology is another interesting topic that is beyond the 
scope of this work (For a most recent study, see Ref. 92]). Besides, as a tentative study 


on the phenomenology of the solar SIDM, our discussion can be further extended to other 
SIDM models, not confined to the case of scalar mediator or Majorana DM. 

Simple as it is, this work has captured the main features of the solar SIDM. To put it in 
a nutshell, the interplay between the long-range dark force and the solar neutrino detection 
is encoded in the following expression up to a possible correction due to the self-capture 
saturation, 


r,ip - 2 


C n + — 


C s 2 K ff 


^ (cTUrel) swave + ( cr ' y rel) pwave 


BR, 


■signal 5 


(41) 


where (crn re i) swave ((0Wei) pwave ) is the s-wave (p-wave) thermal cross section, and BR signal 
is the branch ratio for the signal annihilation channel. Qualitatively speaking, a model 
with a combination of enhanced C s and velocity-suppressed (<ru re i) pwave may maximize the 
novel effects of the SIDM, as long as such enhancement has not been offset by an extremely 
small BR S ig na i. For instance, in the model where 0 decays to the active neutrinos before the 


BBN 66], BR signa i will be significantly boosted if such decay proceeds before 0 reaches the 


Earth. In this case, even for a smaller z the neutrino detector SK will lead in the sensitivity 
on C n . 
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Appendix A: Distribution of the self-interacting dark matter in the sun 


The distribution of the DM particles trapped in the Sun has been well studied in the 
literature [odl 951. While in Ref. 931 authors adopts a Maxwellian velocity distribution 


f x (v) oc exp (—E-k/T x ) = exp [— m x v 2 / (2 T x )],* to describe the non-thermal equilibrium 
between the DM particles and solar nuclei, the studies in Ref. [94[ 95 indicate a deviation 
from the Maxwellian form, in a manner that the actual velocity distribution is suppressed 
at the tail and tends to be anisotropic at large radius. Such deviation can be attributed 
to the fact that the energetic collisions that send the DM particles into the high orbits 
occur predominantly near the hotter core of the Sun. Thus one expects a lower angular 
momentum distribution, or equivalently a larger radial-velocity component (y r /v) for the 
high-energy orbits. This is illustrated in Fig. 2 of Ref. 95], where the normalized average 
kinetic energy of the DM particles (2/3) ( E k) /T x decreases toward the outer layers of the 
Sun, and the average square of the radial velocity component {]( v r /v ) 2 ) deviates from 1/3 at 
radii larger than around 0.2 R Q , indicating an anisotropic nature of the velocity distribution 
beyond that radius. Nevertheless, it should be stressed that while such correction to the 
thermal distribution is necessary for the evaluation of the energy transport and evaporation 
in the Sun that rely sensitively on the high end of the velocity distribution 9J, [95], the 
phenomenology of interest in this study such as self-capture and annihilation can still be 
described with the bulk of the Boltzmann form to a good approximation. 

To illustrate this point, we first investigate the solar DM distribution by following the 
“Brownian motion” method pioneered by the author of Ref. jo3|. A total number of 10 5 
successive collisions between the 10 GeV trial DM particle and solar elements are simulated. 
In the left panel of Fig. [9] we present the quantity (2/3) (l?k) /T x and ((v r /v) 2 ) from the 
solar center up to 0.1 Rq, a radius sufficient to envelope almost all the DM particles trapped 
in the Sun (see Eq. (fT5l) ). The curves would be horizontal (at 1 and 1/3 respectively) for a 
thermal distribution. We compare the simulated velocity distribution with the Maxwellian 
form in the right panel, where the velocity variable v is nondimensionalized in terms of the 
velocity unit (GM q /Rq) 1 ^ 2 » 436 km • s -1 . It is evident that the Maxwellian distribution 
with the characteristic temperature T x indeed provides a good description of the realistic 
solar DM distribution for the mass of interest. T x = 0.967T 0 (O) is determined from Eq. (fT3l) . 
For larger DM masses, one can approximate T x as the solar center temperature T 0 (0). 

In addition, the Boltzmann form is expected to remain an accurate description for the 
solar SIDM because it naturally describes the thermal distribution of the SIDM particles 
gathering in the solar core, while the effective temperature T x accounts for the non-thermal 
equilibrium between the SIDM particles and the solar nuclei to a good approximation. 
Therefore we still adopt the Boltzmann distribution in the evaluation of self-capture and 
annihilation in the Sun. 


* By demand of satisfying the Boltzmann equation, the Maxwellian velocity distribution in turn leads to the 
Boltzmann distribution f x (r, v) oc exp [— (V (r) + EG) /T x \, and hence the density distribution n x (r) oc 
exp [—V (r) /T x ], with V (r) the gravitational potential of the DM particle at radius r. 
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Figure 9: Left panel: the normalized average kinetic energy (2/3) (E^) /T x (solid) and the average 
square of the radial velocity component ^(ty/u) 2 ^ (dotted) for DM mass rn x = lOGeV. Right 
panel: The integrated probability density function (P.D.F) of the solar DM particles obtained from 
the Brownian motion simulation (red histogram). For comparison also shown is the Maxwellian 
distribution (blue lines). See text for details. 


Appendix B: Self-evaporation rate of the trapped dark matter 


In this appendix we will verify the assumption made below Eq. (J5J) that the evaporation 
rate due to the self-interaction of dark matter is negligible compared to the annihilation 
rate. To this end we try to estimate the dark matter self-evaporation rate in the approach 


proposed in Ref. 68], by which we start with the self-evaporation rate e s in a unit volume 


nr 


e s = 


d 3 w 1 f x ( wj) / d 3 w 2 f x (w 2 ) |wi-w 2 | 


d<7 (w 1; W 2 ; Vi, V 2 ) 3 3 

--—-- d°v 1 d c ‘v 2 , 

av\dv2 


(Bl) 


where wq and w 2 are respectively the initial velocities of the trapped DM before scatter¬ 
ing, with their corresponding Maxwellian distributions f x (w x ) and f x (w 2 ) introduced in 
Eq. (l30lh Vi and v 2 being the final velocity after collision, one of which must exceed the 
local escape velocity v esc in the Sun so as to be evaporated, and da (wi, w 2 ; Vi, v 2 ) fdv\dv 2 
schematically represents relevant differential cross section . Then the conservation of energy 
and momentum implies 


/x( W l)/x( W 2) = fx ( V 0 fx ( V 2) , 
|w x - W 2 | = | V i - V 2 | , 


(B2) 


so we have 


nr 


= ~2 rf3 D/x( v i) / d 3 v 2 / x (v 2 ) |vr - v 2 | 


da(y i,v 2 ;wi,w 2 ) 3 3 


dw\dw 2 


-d 6 wid 6 w 2 


~ n 


= n 


d 3v fx ( v ) v 


do (v, 0; w 1; w 2 ) 3 3 


dw\dw 2 


-d w\d w 2 


d 3 v f x (v) wr !c («), 


(B3) 


' v>v e , 
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where we make approximation |vi — V2I ~ |vi| = V\ for the case V\ > v esc u 2 , in which the 
contribution to evaporation rate doesn’t suffer the Maxwellian tail suppression of / x (v 2), 
and cr sc (v) is no other than the self-capture cross section with incident velocity v. On the 
other hand, the annihilation rate in a unit volume is written as 

n\ (cr ann u) = n 2 x j d 3 Vi f x (vi) J d 3 v 2 f x (v 2 ) |v : - v 2 | cr ann (|vi - v 2 |) 

~ n\ J d 3 (y/yfi) f x (v/v^) va ann (v) 

= n\ J d 3 vf x (v) V2 v a ann ( 'y/2v ) , (B4) 

where cr ann (r>) is the annihilation cross section with relative velocity v. Given that compared 
with Eq. (1B4I) the integral over subspace f v>v d 3 v f x (v) in Eq. (1B3D so overwhelmingly 
suppresses e s in the DM mass region of our interest that we can safely ignore the contribution 
from dark matter self-interaction in this study. Finally it is worth mentioning that while our 
discussion is based on the assumption of a large mean free path for the SIDM particles, the 
self-evaporation will be heavily suppressed if the optical depth is much smaller than the DM 
radius. Because in such case, the would-be evaporation events will end up being blocked by 
multiple collisions, which further favors our conclusion. 
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